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ABSTRACT 

We present photometry of the giant extrasolar planet WASP-4b at 3.6 and 4.5 /im taken with the 
Infrared Array Camera on board the Spitzer Space Telescope as part of Spitzer's extended warm 
mission. We find secondary eclipse depths of 0.319%±0.031% and 0.343%±0.027% for the 3.6 and 
4.5 /im bands, respectively and show model emission spectra and pressure-temperature profiles for 
the planetary atmosphere. These eclipse depths are well fit by model emission spectra with water 
and other molecules in absorption, similar to those used for TrES-3 and HD 189733b. Depending on 
our choice of model, these results indicate that this planet has either a weak dayside temperature 
inversion or no inversion at all. The absence of a strong thermal inversion on this highly irradiated 
planet is contrary to the idea that highly irradiated planets are expected to have inversions, perhaps 
due the presence of an unknown absorber in the upper atmosphere. This result might be explained by 
the modestly enhanced activity level of WASP-4b's G7V host star, which could increase the amount 
of UV flux received by the planet, therefore reducing the abundance of the unknown stratospheric 
absorber in the planetary atmosphere as suggested in Knutson et al. (2010). We also find no evidence 
for an offset in the timing of the secondary eclipse and place a 2a upper limit on ecoscj of 0.0024, 
which constrains the range of tidal heating models that could explain this planet's inflated radius. 
Subject headings: eclipses - planetary systems - stars: individual: (WASP-4b) — techniques: photo- 
metric 



1. INTRODUCTION 

Observations of the emergent spectra from transiting 
extrasolar planets with the Spitzer Space Telescope have 
enabled us to probe the atmospheres of a class of giant 
extrasolar planets known as "hot Jupiters" . These plan- 
ets have masses and radii similar to the gas giants in our 
solar system, but orbit very close to their parent stars, 
with equilibrium temperatures ranging from 1000-2500 
K. By measuring the wavelength-dependent decrease in 
light when the planet moves behind the star in an event 
known as a secondary eclipse, we ca n construct a dayside 
emission spectrum for th e planet (|Deming et al.l 120051 : 
ICharbonneau et aT1l2005l) . During its cryogenic mission, 
Spitzer obtained multi-wavelength observations for fif- 
teen extrasolar planets during secondary eclipse. The 
results of these studies indicate that hot Jupiter atmo- 
spheres can be distinguished by the presence or absence 
of a strong temperature inversion in th e upper atmo- 
sphe re (e.g.. iBurrows et al.l I2007U 120081: iFortnev et al.l 
[20081: lBa7manll2008t iMadhusudhan et al.ll2009D . ~ 
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The Spitzer Space Telescope is continuing to survey 
hot Jupiter emission spectra during its post-cryogenic 
mission. After its cryogen was exhausted in May 
2009, only the 3.6 and 4.5 /im channels of the In- 
frared Array Camera (IRAC; Fazio et al. 2004) in- 
strument are operational. Fortunately, these two wave- 
bands are well placed to constrain the range of pos- 
sible models for these atmospheres. Planets with- 
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2009). TrES-1 (ICharbonneau et al.l I2005D and TrES-3 
(jFressin et al.ll2O10t) . are best described by models that 
exhibit H2O and CO absorption features, which cause a 
decrease in the eclipse depth at 4.5 /im relative to 3.6 
/tm. A strong thermal inversion changes these features 
from absorption to emission, therefore increasing the flux 
at wavelengths greater than 4 /im in the atmospheres 
of planets, such as HP 209458b (IDeming et al.l I2005t 
Richardson et al.|[2007l: IKnutson et al.H2008l: ISwain etldl 



2009f ). Of the systems already observed with Spitzer, 
eleven have been found to possess strong temperature 
inversions (see Knutson et al. 2010 for a review). 

In this paper, we present measurements of the tran- 
siting extrasolar planet WASP-4b spanning two times of 
secondary eclipse. WASP-4b is a 1. 24 Mj uv planet or- 
biting at 0.023 AU from a G7V star (jWilson et al.l l2008l: 
IGillon et al.ll200l IWinn et~aT1l2009l . If we assume that 
the planet absorbs all incident flux and re-emits that 
flux as a blackbody from the dayside alone, we calcu- 
late a maximum dayside effective temperature of about 
2000 K. This highly irradiated planet provides an excel- 
lent test case for the correlation between temperature 
inversions and stellar irradiation (for a recent review see 
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Fig. 1. — Background estimate vs. time for 3.6 /im images. We 
estimate the background by fitting a Gaussian to the central region 
of the histogram of counts in the entire array. The background 
estimates exhibit a ramp-like behavior, while also varying between 
three distinct levels. 

Wheatley et al. 2010). It has been hypothesized that 
absorbers such as gas-phase TiO in the upper atmo- 
sphe re trap stellar irradi ation, creating a thermal inver- 
sion ([Hubeny et al .1120031 ) . However, both because TiO is 
a heavy molecule and because titanium can condense into 
solid grains in night-side and day-side cold traps, signifi- 
cant macroscopic mixing would be required to maintain 
it in the upper atmosphere; it is not clear whether such 
vigorous mixing should be expected in a stably strati- 
fied atmosphere (Spie geTet al.1 12009). This theory also 
fails to explain the presence of a temperature inversion 
in XO-lb's atmosphere, as this planet has a dayside tem- 
perature well below the condensation point for TiO. One 
alternative theory suggests that temperature inversions 
could be explained by absorption o f UV and violet vis ible 
light by sulfur-containing species (jZahnle et al.l l2009). 

WASP-4b has a radius of 1.365±0.021 R Jup 
(W inn et all 120091 ) . which is l arger than predicted 
by models of irradiated planets dBurrows et al.l l2007at 
iFortnev et all 120071 : iGuillotl [2001 . placing it among a 
subset of "bloated planets". One possible explanation 
is that the inflated radius is caused by tidal heating 
due to ongoing orbital circularization. Using formulae 
from Liu et al. (2008), Winn et al. (2009) find that 
an orbital eccentricity between 0.002 and 0.02 would 
produce enough heat to inflate the planet to its observed 
size. Using radial velocity measurements, Mahusudhan 
& Winn (2009) find a 95.4% confidence upper limit on e 
of 0.096. By measuring the time of secondary eclipse, we 
can place a much tighter upper limit on the parameter 
ecosw, which will help determine whether tidal heating 
is a viable explanation. 

In Sec. 2 we describe the observations and outline our 
fits to the data. In Sec. 3, we compare our results to the 
predictions of atmospheric models. Finally, in Sec. 4, we 
present our conclusions. 

2. OBSERVATIONS AND METHODS 

We observed a secondary eclipse of WASP-4b in the 
4.5 /xm band on UT 2009 December 6 using IRAC on 
board the Spitzer Space Telescope. We observed in full 
array mode with a 10.4 s integration time, yielding a total 
of 2115 images over a period of 7.7 hr. We observed 
a second secondary eclipse in the 3.6 /jm band on UT 
2009 December 9 using the same 10.4 s integration time, 
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Fig. 2. — Photometry at 3.6 and 4.5 /im vs. time from center 
of secondary eclipse. The decorrelation functions to correct for 
intrapixel sensitivity are overplotted. We use a linear function of 
x and y position to correct the 4.5 fim photometry and a linear 
function of x, y and time for the 3.6 /an observations. 

acquiring 2115 images in 7.7 hr. 

We perform photometry on the Basic Calibrated Data 
(BCD) files produced by version S18.13.0 of the Spitzer 
pipeline. These data files are dark-subtracted, linearized, 
flat fielded and flux-calibrated. The cBCD images have 
been further corrected for artifacts due to bright sources, 
such as column pulldown, but these corrections have an 
unknown effect on time series photometry and we there- 
fore elect to use the standard BCD images in our anal- 
ysis. We extract the UTC-based Julian date for each 
image from the FITS header (keyword DATE.OBS) and 
correct to mid-exposure. We convert to UTC-based BJD 
using the JPL Horizons ephemeris to estimate Spitzer 's 
position during the observations. 

We correct for transient "hot pixels" in a 20x20 pixel 
box around the star by comparing each pixel's intensity 
to the median of the 10 preceding and 10 following frames 
at that position. If a pixel in an individual frame has 
an intensity > 3<r from the median value, its value is 
replaced by the median. We corrected 0.32% and 0.35% 
of the pixels in the box in the 3.6 /im and the 4.5 /im 
band images, respectively. 

We estimate the background by fitting a Gaussian to 
the central region of the histogram of counts in the entire 
array. We find that the background varies significantly 
from frame to frame for both channels. The background 
values, which are plotted in Figure 1 for the 3.6 /im band 
images, display a ramp-like behavior, while also varying 
between three distinct levels. We find a similar pattern in 
channel 2. This behavior is likely a ubiquitous feature of 
warm Spitzer , as it is also observe d in the warm Spitze r 
analysis of CoRoT-1 and CoRoT-2 (|Deming et al.ll2010f) . 

We use three methods to measure the position of the 
star on the array. We calculate the flux-weighted cen- 
troid within 5.0 pixels of the approximate center of the 
star, fit a 2D Gaussian with a fixed width to a 7x7 pixel 
subarray centered on the brightest pix el of the star (e.g., 
lAgol et al. 2010; St evenson et al.ll2010D and fit Gaussians 
to the marginal x and y sums using GCNTRD, which 
is part of the standard IDL astronomy library. Each 
method yields eclipse depths consistent to within la. We 
found that using GCNTRD estimates for channel 1 and 
the 2D Gaussian estimates for channel 2 produced the 
smallest reduced chi-squared for the fits and therefore 
elect to use these position estimate methods. (2D Gaus- 
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Fig. 3. — Photometry at both wavebands after decorrelation vs. 
time from center of secondary eclipse. The data are binned with 
6.6 minute intervals. The error bars are based on the scatter of 
the individual points in each bin. The best-fit eclipse curve is 
overplotted. 

sian fits produced x 2 =2535 and % 2 =2239 for channel 
1 and 2, respectively, whereas GCNTRD fits produced 
X 2 =1973 and x 2 =2273 for channel 1 and 2, respectively.) 
The difference in x 2 produced for the different position 
estimates is very large in channel 1. While we trim ap- 
proximately the same number of points using both meth- 
ods, we find that the GCNTRD positions result in a lower 
level of correlated noise in the final light curve. The rms 
difference between the 2D Gaussian and GCNTRD po- 
sitions are 0.075 pixels in x and 0.229 pixels in y for 
channel 1 and 0.059 pixels in x and 0.145 pixels in y for 
channel 2. These differences are primarily in the form 
of a constant offset; we find that the relative change in 
position calculated using both methods is quite similar. 

We perform aperture photometry with DAOPHOT us- 
ing apertures ranging from 3.0 to 5.0 pixels in half pixel 
intervals. We carried out our fits for each of these aper- 
tures and found that the eclipse depths and times remain 
consistent for apertures between 3.0 and 5.0 pixels. We 
choose an aperture size of 3.5 for our analysis because it 
minimizes both the probability of hot pixels falling within 
the aperture and the root mean square (rms) scatter in 
the data. 

The position of the star varies by 0.50 pixels in x and 
0.46 pixels in y in the 3.6 /im images and by 0.21 pixels 
in x and 0.26 pixels in y in the 4.5 /im images. We dis- 
card any images where the measured flux, x position or y 
position was > 3cr from the median of the twenty frames 
surrounding the image in the time series. We removed a 
total of 10 images (0.47%) and 15 images (0.71%) from 
the 3.6 /xm and 4.5 /im observations, respectively. 

The measured flux from the star varies significantly 
with its po sition on the pixel (e.g. Charbo nneau et all 
2005, 2008). In order to correct for this intrapixel sensi- 
tivity, we fit the data with linear functions of the x and y 
positions. We fit the 4.5 /im data with a linear function 
of the form, 

/ = fo(ci(x - x Q ) + c 2 (y - ?/o ) + c 3 ) (1) 

where / is the flux measured on the array, /o is the orig- 
inal flux of the star, x and y are the positions of the 
star on the array, xq and j/o are the median values of 
x and y over the time series and the constants c\ - C3 
are free parameters. As a check we also try fits to the 



4.5 /im data using a linear function of time instead of 
the x and y variables described above, but this results 
in noticeably poorer fits (x 2 =2429 for the linear fit with 
four d.o.f. (degrees of freedom) and 2239 for the func- 
tion of a; and y including five d.o.f.). We also try a linear 
function of x, y and time, but find the additional time 
term produces only a negligible improvement in the fit 
(x 2 =2235, six d.o.f.). We fit the 3.6 /im data with a 
linear function in x, y, and time. We find that the lin- 
ear fit in time produces a clear improvement in both the 
chi-squared value (% 2 =1973, six d.o.f. and x 2 =2007, five 
d.o.f. for the fits with and without a linear fit for time, 
respectively) and the amount of correlated noise. We 
also try adding quadratic terms in x and y which are 
usually required when the star falls on the peak of the 
intrapixel curve (center of the pixel). However, we find 
that adding additional degrees of freedom in x and y has 
a negligible effect on the final time series, eclipse values, 
and chi-squared (x 2 =1968, eight d.o.f.) and therefore 
elect to use the linear fit. Figure 2 shows the photome- 
try with the decorrelation functions overplotted for each 
waveband. 

We use a Markov Chain Monte Carlo method (jFord I 
l2005tlWinn et al.ll2007t ) with 10 6 steps to simultaneously 
determine the transit depth, timing of the eclipse, and 
the corrections for intrapixel sensitivity. We use five free 
parameters in the 4.5 /im data and six free parameters in 
the 3.6 /im data, including the linear term in time. We 
set the system parameters (planetary and stellar radii, 
orbital period, and orbital inclination) to the values given 
in Winn et al. (2009). We calculate the eclipse curve us- 
ing the equations from Mandel & Agol (2002). The un- 
certainty for each point is set equal to the rms deviation 
of the out-of-eclipse data after removing the intrapixel 
effect. We also trim the first half hour of data from both 
the 3.6 and 4.5 /im time series because it exhibits larger 
deviations in position, perhaps due to settling of the tele- 
scope at a new pointing. 

We take the median value of the distribution for each 
parameter as our best-fit solution. We calculate symmet- 
ric error bars about the median by finding the range over 
which the probability distribution contains 68% of the 
points above and below the median. The distributions 
for all parameters are nearly Gaussian and there are no 
strong correlations between parameters. Best-fit eclipse 
depths and times are shown in Table 1. As a check, we 
ran a second independent Markov chain for each chan- 
nel and obtained identical results. Figure 3 shows the 
photometry after it has been corrected with the best-fit 
intrapixel correlation function with the best-fit eclipse 
curves overplotted. 

We als o calculate error bars using the 'prayer-bead' 
method (jGillon et al.ll2009f ). We divide our time series 
by the best-fit solution from the Markov chain, shift the 
time series in one point increments, multiply the best-fit 
solution back in and calculate the eclipse depth and time 
for the new data set. The prayer-bead distributions gave 
error bars that were consistent with the Markov chain 
errors and we elect to use the larger of the two errors in 
each case. For channel 1, we use the prayer-bead error for 
the eclipse depth (0.031%) instead of the Markov error 
(0.019%), whereas we use the Markov error for the time 
(1.3 min) instead of the prayer-bead error (0.72 min). 
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TABLE 1 

Summary of Secondary Eclipse Results 



Wavelength (/im) Center of Eclipse (BJD) Depth (%) Eclipse Offset (min) T bright (K) a 

3.6 2455174.87731 ± 0.00087 0.319 ± 0.031 0.5 ± 1.3 1832 ± 71 

4.5 2455172.2011 ±0.0013 0.343 ± 0.027 0.1 ± 1.9 1632 ± 56 



a We calculate the brightness temperature of the planet by finding the flux-weight ed average of the planet-star flux ratio over each Spitzer 
bandpass. We use a 5500 K PHOENIX NextGen model (Hauschildt ct al. 1999) for the stellar spectrum and set the planet's emission 
spectrum equal to a blackbody, then solve for the temperature at which the planet-star flux ratio equals the observed eclipse depth. 



For channel 2, we used the prayer-bead errors for both 
the eclipse depth (0.027%) and time (1.9 min). The cor- 
responding Markov errors are 0.023% and 1.4 min. 

We find that the rms variation in our light curve after 
correcting for intrapixel sensitivity is 1.1 and 1.2 times 
the predicted photon noise from the star at 3.6 and 4.5 
/im, respectively. The reduced chi squared for our fits 
are 1.01 (x 2 =1973, 1945 points, six d.o.f.) and 1.16 (x 2 = 
2239, 1941 points, five d.o.f.) for the 3.6 and 4.5 /mi light 
curves, respectively. The error used in the fits is based on 
the rms deviation of the out-of-eclipse light curve rather 
than the predicted photon noise. Since we use the rms 
error estimates, the reduced \ 2 should theoretically be 
equal to 1.0 if the noise is purely gaussian. The fact 
that we find reduced \ 2 values exceeding 1.0 reflects the 
correlated noise present in our light curves which we take 
into account with the prayer-bead analysis. 

3. DISCUSSION 
3.1. Orbital Eccentricity 

The timing of the secondary eclipse is very sensitive 
to the planet's orbital eccentricity. Assuming a circular 
orbit and accounting for 23.4 second s that light takes 
to travel across the orbit (jLoebl 120051 ) . we would expect 
to see the secondary eclipse occur at a phase of 0.5002. 
In the event that there is significant advection of energy 
to the planet's night side we would expect an additional 
delay due to an offset hot spot on the planet's day side 
causing a change in the shape of ingress and egress. We 
estimate the maximum value of this delay to be 41 sec- 
onds based on a model in which the longitudinal advec- 
tion time is 60% of the radiative time, corresponding to a 
hot region shifted 30 degrees east of the su bstellar point 
([Williams et al J 12005 ICowen fc Agp]|l20Toh . We can use 
the difference between the predicted and observed orbital 
phase of secondary eclipse, including the light travel time 
but neglecting the unknown delay from a nonuniform 
surface brightness, to constrain ecosw, where e is the 
orbit al eccentricity and u is th e argument of pericenter 
(e.g.. lCharbonneau et al.|[2005ft . 

We find that the eclipse is offset from the predicted 
time based on the ephemeris from Winn et al. (2009) 
by 0.5±1.3 and 0.1±1.9 minutes in the 3.6 and 4.5 
/im bands, respectively. We take the average of these 
two values weighted by the inverse of the variance and 
find a mean of 0.4±1.0 min, corresponding to ecosw = 
0.00030 ±0.00086. We place a 2a upper limit on |ecosw| 
of 0.0024, where we have calculated this limit by integrat- 
ing over the histograms for the eclipse time. We integrate 
over the histograms from the Markov chain distribution 
for channel 1 and the prayer-bead distribution for chan- 
nel 2. This upper limit implies that unless our line of 
sight happens to align very closely with the planet's ma- 



jor axis (i.e. the argument of pericenter u is close to tt/2 
or 37r/2) the orbit is nearly circular. 

Ibgui, Burrows, & Spiegel (2010) investigate the extra 
core power that would be needed to explain the otherwise 
anomalously large radius of WASP-4b. They find that 
approximately 7.8 x 10~ 8 L of heating would be nec- 
essary for solar-metallicity opacity atmospheres, which 
decreases to 10 -8 L Q for 10 x solar opacity atmospheres. 
Less power is necessary if the atmosphere helps retain 
more heat, as in the 10 x solar case. If this heating is 
due to tides, and the eccentricity is on the order of 0.001 
and is maintained by an external planetary perturber in 
the system (as yet unidentified; Mardling 2007), then the 
Q' tidal dissipation parameter would be roughly between 
3 x 10 4 and 2 x 10 5 . In the Ibgui, Burrows, & Spiegel 
(2010) paper, a value of 0.096 is assumed for the eccen- 
tricity and this leads them to derive a "best-estimate" 
range for Q' between 3 x 10 8 and 2 x 10 9 . With our 
new constraint on the eccentricity of WASP-4b's orbit, 
and using the calculations of Ibgui, Burrows, & Spiegel 
(2010), we now obtain a range for Q' that is more in 
line with the measured value of Jupiter of 10 5 — 10 6 
(|Goldreich fc Soterj[l96l lYoder fc Pealdll98lh . 

3.2. Atmospheric Temperature Structure 

In this paper we examine two distinct classes of hot 
Jupiter models. Figure 4 shows our planet/star contrast 
ratios and three models for the planetary atmosphere 
derived from one-dimensional, plane-parallel atmosphere 
codes following Fortney et al. (2008). One model as- 
sumes the presence of the absorber TiO in the upper 
atmosphere at equilibrium abundances, whereas the two 
remaining model atmospheres contain no TiO. Fortney 
et al. (2008) parameterize the unknown redistribution of 
energy to the planet's nightside by varying the stellar 
flux incident at the top of the planetary atmosphere by a 
geometric factor to account for dayside average (/=0.5) 
or planet- wide average (/=0.25) conditions. The slope 
between the 3.6 and 4.5 /im points on the model with 
TiO (green) is much too steep to fit both measurements 
simultaneously. Wc find that WASP-4b is best fit by the 
orange model with no TiO (no inversion) and geometric 
factor /=0.60, resulting in a very hot dayside. This is a 
reasonable choice, as the projected area of the substellar 
point (/=1) is maximized during secondary eclipse while 
contributions from the cooler regions near the day-night 
terminator are correspondingly re duced, giving an aver - 
age value of /=2/3 at opposition (jBurrows et al-l feOOS). 
The dayside pressure-temperature profiles for these three 
models are displayed in Figure 5. 

Figure 6 shows three models for the planetary atmo- 
sphere with greater degrees of freedom following Burrows 
et al. (2008). While the Fortney et al. models con- 
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Fig. 4. — Dayside planet/star flux ratio vs. wavelength for three model atmospheres (Fortncy ct al. 2008) with the band-averaged flux 
ratios for each model superposed (squares). The measured contrast ratios are overplotted (black circles). One model (green) represents 
an atmosphere containing TiO in the upper atmosphere at equilibrium abundance. The other two models (orange and magenta) contain 
no TiO. The parameter / represents the redistribution of energy over the planet's surface, where /=0.50 corresponds to dayside only 
redistribution and /=0.25 corresponds to uniform redistribution over the entire planet. We obtain the best fit to our measurements by a 
model with no TiO and little redistribution (orange). 
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Fig. 5. — Dayside pressure-t emperature profiles for the three 
model atmospheres in Figure 4 (Fortncy ct al. 2008). The green 
model contains TiO in the upper atmosphere and exhibits a strong 
temperature inversion for pressures below 0.01 bars. The orange 
and magenta profiles represent atmospheres with no TiO but have 
different values of the redistribution parameter /. The /=0.25 
model has full redistribution of energy to the nightside, resulting 
in a cooler dayside profile. The hotter /=0.60 model provides the 
best fit to our measurements. We also indicate the approximate 
locations of the 3.6 and 4.5 fim photospheres (solid squares) for 
each model, estimated here as the pressure at which the model 
temperature matches the measured brightness temperature in each 
bandpass. 

tain TiO in either zero or equilibrium abundance, the 
Burrows et al. models contain an unknown absorber at 
various optical opacities, parameterized by n e in cm 2 /g. 
Burrows et al. (2008) add a heat sink at a pressure range 
of 0.01 to 0.1 bars to model energy redistribution from 
the day to the nightside. As energy is most likely redis- 
tributed deep in the planetary atmosphere, this method 
for modeling heat transfer is physically motivated, but 
contains more degrees of freedom than the Fortney et 
al. models. General circulation models for these plan- 
ets indicate that redistribution occurs continuously over 



a range of pressures (e.g.. iShowman et al.ll2008l ). and we 
note that the range of pressures selected for our parame- 
terized redistribution model can have a modest effect on 
the resulting pressure-temperature profiles, although it 
does not affect our main conclusions in this paper. The 
dimcnsionless parameter P n is a measure of the day to 
nightside energy redistribution, where P„=0.0 represents 
no redistribution and P n =0.5 represents full redistribu- 
tion to the nightside. 

Burrows et al. use a 55 00 K Kurucz atmospher e 
model for the stellar spectrum (jKuruc 3GI11GI94, 2005), 
wherea s Fortney et al. use a 55 00 K PHOENIX NextGen 
model (|Hauschildt et al]|1999f) . As a check, we calculate 
the Burrows et al. planet-star flux ratio models using a 
PHOENIX NextGen stellar spectrum instead of the Ku- 
rucz spectrum and find that the differences are minimal 
and comparable to the differences caused by the uncer- 
tainty in the star's effective temperature. We find that 
the differences in the band-integrated flux-ratios using 
the two different stellar models vary between 0.003 and 
0.005% and are therefore negligible when compared to 
our measurement errors. 

We show an inverted atmosphere model (red), with n e 
and P n set equal to the best-fit va lues for the archetype 
inver ted atmosphere HD 209458b (jBurrows et al.ll2007bl 
2008). This inverted model is a poor fit to our mea- 
sured contrast ratio at 4.5 fim. We find the best match 
is the model with a small amount of stratospheric ab- 
sorber with K e =0.03 cm 2 /g and relatively efficient day- 
night circulation with P„=0.3. The band-integrated flux 
ratios for this model (green) fall within la of the mea- 
sured ratios in both bands. The pressure-temperature 
profiles in Figure 7 show that this best-fit model exhibits 
a modest temperature inversion for pressures below 0.01 
bars, much weaker than the archetype inverted atmo- 
sphere HD 209458b. The blue non-inverted atmosphere 
model with parameters K e =0.0 and P„=0.1 fails to fit 
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Fig. 6. — Dayside planet/star flux ratio vs. wavelength for three model atmospheres (Burrows ct al. 2008) with the band-averaged 
flux ratios for each model superposed (squares) to account for the widths of the Spitzer bandpasses. The measured contrast ratios 
are overplotted (black circles). The blue model represents a non-inverted atmosphere (k e =0.0) with redistribution parameter P n =0.1. 
An inverted atmosphere model is shown in red (ft e =0.1), which exhibits water features in emission instead of absorption. The green 
model represents an atmosphere with a small amount of upper-atmosphere absorber, with optical opacity K e equal to 0.03. The Spitzer 
measurements are best matched by this model, which suggests that the atmosphere of WASP-4b has a moderate thermal inversion in its 
upper atmosphere. 



our measurements at both wavebands. 

We note that while the best-fit Burrows et al. model 
indicates moderately efficient [P n = 0.3) day-night cir- 
culation, the best-fit Fortney et al. model with / = 0.60 
requires minimal day-night circulation. It is perhaps not 
surprising that these relatively simple models disagree, 
given the differences in their treatment of the incident 
flux, optical opacities, and energy loss (if any) to the 
night side. We find some tentative evidence that this 
disagreement m ay be systematic, as published results 
for HP 189733b (IKnutson et al.l l2009ah and HD 209458b 
(|Fortnev et al.ll2010f ) from Fortnev et al. favor a hot day- 
side whereas Burrows et al. models predict greater en- 
ergy redistribution (P„ = 0.30 and 0.15 for HD 209485b 
and HD 189733b, respectively; Burrows et al. 2007b, 
2008). Multi- wavelength phase curve observations allow 
us to test these predictions, at least for the brightest 
systems (e.g., Knutson et al. 2007, 2009a). 

Spitzer infrared observations indicate that the atmo- 
spheres of extrasolar giant planets tend to exhibit prop- 
erties ranging between two differing types, exemplified 
by HD 189733b, whose emission spectrum features wa- 
ter and other molecules in absorption, and HD 209458b, 
which exhibits these features in emission. Table 2 shows 
the published values of K e and P n for a range of plan- 
ets with Spitzer observations; WASP-4b is similar to HD 
189733b and TrES-3 in that it requires a relatively small 
amount of absorber as compared to HD 209458b. 

Assuming WASP-4b absorbs with zero albedo and re- 
emits on the dayside only, the planet's predicted day- 
side effective temperature is approximately 2000 K. If 
the planet emits uniformly over both hemispheres, we 
would expect an effective temperature of about 1650 K. 
We fit both measured eclipse depths simu ltaneously using 
a 550 K PHOENIX NextGen model (jHauschildt et al.l 
1999) for the stellar spectrum and a blackbody for the 
planet's spectrum, and find that WASP-4b has a best-fit 



2500 




10 2 10° 
Pressure (bar) 

Fig. 7. — Dayside pressure-temperature profiles for three model 
atmospheres with va rious values of the parameters P n and re e 
(Burrows ct al. 2008). The blue model represents an atmosphere 
with no inversion. The red model corresponds to an atmosphere 
with an additional absorber with optical opacity He =0.1 cm 2 /g. 
The absorber, which is added high up in the atmosphere where 
the pressure is below 0.03 bars, traps stellar irradiation and cre- 
ates a temperature inversion. The green model with re e =0.03 and 
P n =0.3 provides the best fit to our measurements of WASP-4b. 
This model exhibits a slight temperature inversion for pressures 
less than 0.01 bars. Burrows et al. (2008) add a heat sink at a 
pressure range of 0.01 to 0.1 bars to model energy redistribution 
from the day to the nightside, which contributes to the decrease 
in dayside temperatures between 0.05 and 1.0 bars for the P n =0.3 
models. We also indicate the approximate locations of the 3.6 
and 4.5 fim photospheres (solid squares) for each model, estimated 
here as the median pressure of the r =2/3 surface over the range of 
wavelengths spanned by each bandpass. We find the same approx- 
imate photosphere locations by solving for the pressure at which 
the the temperature of the model matches the measured bright- 
ness temperature in each band. Due to the width of the Spitzer 
bandpasses, we actually see flux from a wide range of pressures. 
Typical ranges are 7xl0~ 3 - 2xl0 _1 bars and 2xl0~ 4 - lxlO -1 
bars at 3.6 and 4.5 (im, respectively. 
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TABLE 2 

Effective Temperature and Burrows et al. Model 
Parameters for Extrasolar Giant Planets 



Name 


Teff (K) a 






Reference 


TrES-3 


2000 


0.01 


0.3 


Fressin et al. 2010 


WASP-4b 


2000 


0.03 


0.3 


this paper 


HD 189733b 


1400 


0.035 


0.15 


Grillmair et al. 2008 


HD 209458b 


1700 


0.1 


0.3 


Burrows et al. 2007b 


TrES-4 


2100 


0.1 


0.3 


Knutson et al. 2009b 


XO-lb 


1400 


0.1 


0.3 


Machalek et al. 2008 


XO-2b 


1600 


0.1 


0.3 


Machalek et al. 2009 


TrES-2 


1800 


0.3 


0.3 


Spiegel & Burrows 2010 


HAT-P-7b 


2500 


1.1 


0.0 


Spiegel & Burrows 2010 



jor axis, the planet's orbit must be nearly circular. Al- 
though this upper limit does not rule out tidal heating, 
it constrains the range of tidal heating models that could 
explain this planet's inflated radius. 

We find secondary eclipse depths of 0.319%±0.03f % 
and 0.343%±0.027% for the 3.6 and 4.5 ^m bands, re- 
spectively. These results are consistent with a spectrum 
exhibiting water and CO in absorption. We find that the 
atmosphere can be well characterized by models with a 
modest or no thermal inversion. Measurements at other 
wavelengths would help to distinguish between these two 
models. The absence of a strong thermal inversion makes 
WASP-4b an exception to the rule that inversions are 
found on planets that receive higher stellar irradiation. 
Other exceptio ns include the highl y irradiated extrasolar 
planet TrES-3 (Fressin et al. 201(3) which does not have a 
temperature inversion and XO-lb (jMachalek et al.l[2008| ) 
which possesses a temperature inversion despite being 
relatively cool. These planets indicate that there must 
exist additional stellar or planetary parameters, other 
than equilibrium temperature, responsible for determin- 
ing the relative strengths of thermal inversions in hot 
Jupiter atmospheres. 

This work demonstrates that Warm Spitzer, which 
operates with the 3.6 and 4.5 ^m channels only, can 
be successfully used to characterize the properties of 
hot Jupiter atmospheres. The increasing availability 
of ground-based ec l ipse detections in t he near-IR (e.g., 
Gillo n et al.1 [200l iCroll et all l2010al|bl : iGibson et all 
20101 : iLopez-Morales et al.ll20 10) will also help to resolve 
ambiguities in the interpretation of the Spitzer data for 
many of these planets. Indeed, our models predict that 
WASP-4b should have an eclipse depth of 0.1-0.2% in 
the Ks band (2.15 microns). Croll et al. (2010b) mea- 
sured a secondary eclipse depth of 0.1331QQU in this 
same bandpass for TrES-3, which has an apparent bright- 
ness and other properties similar to those of WASP-4. 
The TrES-3 K-band detection augers well for a similar 
WASP-4b measurement, which would provide a further 
point of comparison for the atmospheric models that we 
present. By the end of its post-cryogenic mission, Spitzer 
will observe more than twenty systems during secondary 
eclipse. When combined with the nineteen systems ob- 
served during the cryogenic mission, as well as any avail- 
able ground-based detections, these results will allow us 
to search for correlations with other system parameters 
that could provide valuable clues to the origin of temper- 
ature inversions in hot Jupiter atmospheres. 

This work is based on observations made with the 
Spitzer Space Telescope, which is operated by the Jet 
Propulsion Laboratory, California Institute of Technol- 
ogy, under contract with NASA. Support for this work 
was provided by NASA. Heather A. Knutson is sup- 
ported by a fellowship from the Miller Institute for Basic 
Research Science. Eric Agol acknowledges the support 
of NSF CAREER Grant No. 0645416. 



a Predicted blackbody temperature for the planet assuming an 
albedo of zero and no nightside redistribution of energy. 

blackbody temperature of 1700 K. Given such high irra- 
diation, it is somewhat surprising that WASP-4b exhibits 
at most a relatively weak thermal inversion. WASP-4b 
is therefore an exception to the general trend that highly 
irradiated planets are more likely to have strong thermal 
inversions. 

In Knutson et al. (2010) we propose that there ex- 
ists a correlation between temperature inversions and 
the activity levels of the host star, where the increased 
UV flux from active host stars destroy the compounds 
that are responsible for producing temperature inver- 
sions. We use Ca II H & K line strengths as indicators 
of stellar activity levels. In Knutson et al. (2010) we 
obtain Keck HIRES spectra for WASP-4b and find the 
Ca II H & K line strength estimates are Sjyif =0.194 and 
log {Rhk 1 )—— 4.865, assuming a model B — V color of 
0.74 for a 5500 K star. These line strengths indicate that 
WASP-4 is a moderately active star, with a log {Rhk 1 ) 
value that falls near the division between classes. How- 
ever, WASP-4b's smaller orbital distance relative to HD 
189733b means that it intercepts proportionally more of 
its star's flux, and as a result we estimate that the UV 
flux per unit area incident at the surface of WASP-4b 
is approximately half that received by the planet HD 
189733b and twice that received by WASP-2 (see discus- 
sion in Knutson et al. 2010). 

We also calculate a value for the empirical index de- 
fined in Knutson et al. (2010) as the difference be- 
tween the slope across the measured 3.6 and 4.5 /im 
eclipse depths and the slope of the best-fit blackbody 
function for the planet, which provides an observational 
means to distinguish between the two hot Jupiter at- 
mosphere types. We find a value of — 0.09±0.04 in this 
index for WASP-4b, which suggests that this planet is 
best classified in the same type as HD 189733b (index of 
-0.15±0.02) and TrES-3b (index of -0.10±0.05). Plan- 
ets with strong inversions typically have positive values 
in this index, therefore this result is consistent with our 
earlier conclusion that WASP-4 displays at most a rela- 
tively weak temperature inversion. 

4. CONCLUSIONS 

We observed two secondary eclipses of the extrasolar 
planet WASP-4b at 3.6 and 4.5 /im as part of Spitzer's 
extended warm mission. By measuring the time of the 
eclipse, we estimate a 2er upper limit on the parameter 
|ecosw| of 0.0024. This limit implies that unless our 
line of sight happens to align closely to the planet's ma- 



8 



Beerer et al. 



REFERENCES 



Agol, E., et al. 2010, ApJ submitted 
Barman, T. S., 2008, ApJ, 676, L61 

Burrows, A., Hubeny, I., Budaj, J., Hubbard, W. B. 2007a, ApJ, 
661, 502 

Burrows, A., et al. 2007b, ApJ, 668, L171 
Burrows, A., Budaj, J., Hubeny, I. 2008, ApJ, 678, 1463 
Charbonneau, D., et al. 2005, ApJ, 626, 523 
Charbonneau, D., et al. 2008, ApJ, 686, 1341 
Cowen, N. B., & Agol, E. 2010, ApJ submitted 
Croll, B., et al. 2010a, ApJ, 717, 1084 
Croll, B., et al. 2010b, ApJ, 718, 920 
Deming, D., et al. 2005, Nature, 434, 740 

Deming, D., Harrington, J., Seager, S., & Richardson, L. J. 2006, 

ApJ, 644, 560 
Deming, D., et al. 2010, Nature, in prep 
Fazio, G. G., et al. 2004, ApJ, 154, 10 
Ford, E. 2005, ApJ, 129, 1706 

Fortney J. J., Marley, M. S., Barnes, J. W. 2007, ApJ, 659, 1661 
Fortney J. J., Lodders, K., Marley, M. S., & Freedman, R. S. 

2008, ApJ, 678, 1419 
Fortney J. J., et al. 2010, ApJ, 709, 1396 
Fressin, F., et al. 2010, ApJ, 711, 374 
Gibson, N. P., et al. 2010, MNRAS, 404, L114 
Gillon, M., et al. 2009, A&A, 496, 259 
Gillon, M., et al. 2010, A&A, 506, 359 
Goldreich, P., & Soter, S. 1966, Icarus, 5, 375 
Grillmair, C. J., et al. 2007, ApJ, 658, L115 
Grillmair, C. J., et al. 2008, Nature, 456, 767 
Guillot, T. 2008, Physica Scripta Volume T, 130, 014023 
Hauschildt, P. H., et al. 1999, ApJ, 525, 871 

Hubeny, I., Burrows, A., Sudarsky, D., et al. 2003, ApJ, 594, 1011 
Ibgui, L., Burrows, A., & Spiegel, D.S. 2010, ApJ, 713, 751 
Knutson, H. A., et al. 2007, Nature, 447, 183 
Knutson, H. A., et al. 2008, ApJ, 673, 526 



Knutson, H. A., et al. 2009a, ApJ, 690, 822 

Knutson, H. A., et al. 2009b, ApJ, 691, 866 

Knutson, H. A., Howard, A. W., & Isaacson, H. 2010, in prep 

Kurucz, R. 1979, ApJS, 40, 1 

Kurucz, R. 1994, Solar Abundance Model Atmospheres for 0, 1, 

2, 4i an d 8 kra/s CD-ROM 19 (Smithsonian Astrophysical 

Observatory, Cambridge, MA, 1994) 
Kurucz, R. 2005, in Memorie Delia Societa Astronomica Italiana 

Supplement, v. 8, p. 14 
Liu, X., Burrows, A., & Ibgui, L. 2008, ApJ, 678, 1396 
Loeb, A. 2005, ApJ, 623, L45 
Lopez-Morales, M., et al. 2010, ApJ, 716, L36 
Machalek, P., et al. 2008, ApJ, 684, 1427 
Machalek, P., et al. 2009, ApJ, 701, 514 
Mahusudhan, N., & Seager, S. 2009, ApJ, 707, 24 
Mahusudhan, N., & Winn, J. N. 2009, ApJ, 693, 784 
Mandel, K., & Agol, E. 2002, ApJ, 580, L171 
Mardling, R.A. 2007, MNRAS, 382, 1768 
Richardson, L. J., et al. 2007, Nature, 445, 892 
Showman, A. P., et al. 2008, ApJ, 682, 559 
Spiegel, D. S., Silverio, K., & Burrows, A. 2009, ApJ, 699, 1487 
Spiegel, D. S., & Burrows, A. 2010, in prep 
Stevenson, K. B., et al. 2010, Nature, 464, 1161 
Swain, M. R., et al. 2009, ApJ, 690, L114 
Wheatley, P. J., et al. 2010 in prep 
Williams, P. K. G., et al. 2006, ApJ, 649, 1020 
Wilson, D. M., et al. 2008, ApJ, 675, L113-L116 
Winn, J. N., Holman, M. J., & Fuentes, C. I. 2007, AJ, 133, 11 
Winn, J. N., Holman M. J., Carter J. A., et al. 2009, AJ, 137, 

3826 

Yoder, C. F., & Peale, S. J. 1981, Icarus, 47,1 
Zahnle, K., et al. 2009, ApJ, 701, L20 



